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We study thermodynamics of cosmological models in the Horava-Lifshitz theory of gravity, and 
systematically investigate the evolution of the universe filled with a perfect fluid that has the equation 
of state p = wp, where p and p denote, respectively, the pressure and energy density of the fluid, 
and w is an arbitrary real constant. Depending on specific values of the free parameters involved 
in the models, we classify all of them into various cases. In each case the main properties of the 
evolution are studied in detail, including the periods of deceleration and/or acceleration, and the 
existence of big bang, big crunch, and big rip singularities. We pay particular attention on models 
that may give rise to a bouncing universe. 
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I. INTRODUCTION 

Recently, Horava proposed a very attractive quantum 
gravity theory [ij, motivated by the Lifshitz theory in 
solid state physics [2l|, for which the theory is usually 
referred to as the Horava-Lifshitz (HL) theory. In the 
IR limit, the standard general relativity is recovered. By 
construction, it is non-relativistic and UV-renormalizable 
at least around the flat space. The effective speed of 
light diverges in the UV, and may potentially resolve the 
horizon problem without invoking inflationary scenario. 
Since the theory is brand new, detailed studies are highly 
demanded, before any definitive conclusions are reached, 
although a great deal of efforts have already been de- 
voted to these subjects, including the studies of cosmol- 
ogy 0, 3 , and black hole physics , among others Q . 
In particular, in 0,0 the general field equations were de- 
rived. When applying them to cosmology, the complete 
set of field equations were given explicitly, from which 
it can be seen that the spatial curvature is enhanced by 
higher-order curvature terms, and this may allow us to 
address the flatness problem, and provide a bouncing cos- 
mology [H, [lo| • It was also shown that almost-invariant 
super-horizon curvature perturbations can be produced 



sup e 

0. 



However, despite these attractive features, the the- 
ory has already been facing some challenging questions 
[13, [3 ■ In particular, it was shown that the HL theory 
may suffer strong coupling problems due to the break- 
ing of diffeomorphism invariance [l^ . As pointed out in 
14l |. these problems might be solved by preserving the 
prejactablity condition, as was done originally by Horava 
On the other hand, there are a couple of reasons to 
abandon the detailed balance condition. One is due to 
the fact that matter is not UV stable with this condition 
0. It also requires a non-zero (negative) cosmological 
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constant in order to have a correct coupl ing, and breaks 
parity in the purely gravitational sector [l5|. As shown 
explicitly in [l5| , the general theory can be properly for- 
mulated without the "detailed balance" condition, but 
still keeping the projectability condition and preserving 
parity. 

In this paper, we shall focus on the thermodynamics of 
cosmological models in the HL theory of gravity without 
detailed balance, and systematically investigate cosmo- 
logical models for a perfect fluid with the equation of 
state p = wp, where p and p denote, respectively, the 
energy density and pressure of the fluid, and w is an ar- 
bitrary real constant. We shall classify all these models 
according to the values of the parameters involved in the 
models, and study the evolution of the universe for each 
model. By doing so, we shall study the spacetime singu- 
larities, such as the big bang, big crunch and big rip, and 
identify the period(s) when the universe is accelerating or 
decelerating. We pay particular attention on models that 
give rise to a bouncing universe. Speciflcally, the paper 
is organized as follows: In Sec. II, we give a brief intro- 
duction to the HL theory, while in Sec. Ill, we present 
the Friedmann-like field equations. In Sec. IV, we study 
thermodynamics of cosmological models, and in Sec. V. 
we investigate the evolution of the universe when filled 
with a perfect fluid with the equation of state p = wp. 
We study these solutions case by case, and deduce the 
main properties of each model of the universe. Finally, 
in Sec. VI, we present our main conclusions. 

It should be noted that classiflcation of a (non- 
relativistic) matter coupled with a dark energy was con- 
sidered recently in [l6l |. in the framework of Einstein's 
theory, and the corresponding Penrose diagrams were 
presented. Similar considerations were also carried out 
in a series of papers, and particular attention was paid 
to obtain an effective potential V{a) by fltting observa- 
tional data sets [1'^]. In [l8j, such studies were general- 
ized to a perfect fluid with the equation of state p — wp. 
In this paper, we shall generalize these studies to the HL 
cosmology. 
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II. THE HORAVA-LIFSHITZ GRAVITY 
THEORY 

In this section, we shall give a very brief introduction 
to the HL Theory. For detail, we refer readers to P, 0, 
The dynamical variables are N, Ni and gij {i, j = 
1, 2, 3), in terms of which the metric takes the ADM form, 

ds^ = -N'^dt^ + grj {dx' + N'dt) {dx^ + N^dt) , (2.1) 

where TV* = g^^Nj, and the coordinates (i,a;') scale as, 

t £h, x' ix\ (2.2) 

Under the above scaling, the dynamical variables scale as 

N-^N, ^ .9,,, N, fN„ ^ i-^N\ (2.3) 

The total action of the HL theory consists of three 
parts, the kinetic part, 5*^, the potential part, Sy, and 
the matter part, 5*^, given by 

Stotal = + Sy + Sjn 

dtdx'Ny^{Ck + Cy + £„0, (2.4) 



where g is the determinant of the three- metric gij , Cm = 
Cjn (N, Ni, gij, $) the Lagrangian density of matter 
fields, denoted collectively by $, and 

Ck = a {K,,K'^ - XK^) , 



V9 



(2.5) 



where e*-''^ is the antisymmetric tensor with e^^^ = 1, Vfe 
denotes the covariant derivative with respect to gij. Rij 
is the Ricci tensor of the three-metric gij, R = g^^ Rij, 
and Cij and Kij are, respectively, the Cotton tensor and 
extrinsic curvature, defined by 



^ —V, [Ri -6iR 
K^J ^ ^ (ft, - V.iV, - V,iV,) , (2.6) 

where ijij = dgij/dt. The constants a. A, /3,7, C, rj, ^ 
and a are coupling constants. Under the "detailed- 
balance" conditions, they arc not independent, and are 
given by 



2 k2 . 

AtV2(l~4A) ^ _ K^^A 
~ 32(1 -3A) ' 8(1 -3A)' 
_ _ 3^VA^ 
^ "8(1-3A)' 



(2.7) 



where [= SttG/c^] and A are, the Einstein coupling 
and cosmological constants, respectively, and A, /i and 
Lo are the three independent coupling constants of the 
theory. As pointed out in [17,] , one can make an analytical 
continuation of the parameters /i and by 



fi ifi, a; — > —luj , 
so that the coupling constants change as, 

a a, ^ -/?, 7 ^ -7, C - 



-c, 



(2.8) 



(2.9) 



In this paper, we shall not impose the "detailed- 
balance" conditions given by Ea. (|2.7p . so that all the 
constants appearing in the Lagrangian densities given by 
Eq. p.Sp are independent and otherwise arbitrary, subject 
to the constraint. 



a(3A - 1) > 0, 



(2.10) 



a condition that will be clear when we study cosmological 
models in the next section. 

In the IR limit, all the quadratic terms of Rij are 
dropped out, and the total action reduced to 



Stotal ^ J dtdx'N^ [a {K^jK'^ - XK^ 
+rjR^ + ^R + (7], 
which will reduce to the Einstein-Hilbert action, 
1 



Seh — 



IGttG 
by setting x^ — ct, 

X 

IGttG 



■(i?4[ff] 



2A 



EH 



(2.11) 



(2.12) 



A 



a ' 

EH 



where 



500 

9i3 



-~N^ + g'^N,N,, go. 



N,. 



(2.14) 



Note that Condition (|2.10p . together with the one that 
c is real, requires A < 0. To get a positive A, one can 
invoke the analytical continuation of the parameters fi 
and w^, given by Eqs.JlJl) and ([2^]) . 



III. COSMOLOGICAL MODELS IN THE 
HORAVA-LIFSHITZ THEORY 

The homogeneous and isotropic universe is described 
by the metric, 



dr^ 



1 — fcr^ 



(3.1) 
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where cPQ = dO'^+sin^ 6d(j)'^, and k = 0,±1. For a perfect 
fluid, 



Tab = (P+P) UaUb + pgab, 



(3.2) 



where Ua = denotes the four- velocity of the fluid, 
the field equations of the Horava-Lifshitz theory can be 
costed in the forms, 



3a(3A- ^ p-cr 



6k^ 



12fc2 (( + 2ri) 



p + 3H{p+p) = 0, 

where H — a/a. 

Introducing the following quantities, 



(3.3) 
(3.4) 



SttG = 



a(3A - 1) 



, PA = -PA = -0-, 



Pk = -3pk Pdr^ 3pdr = (3.5) 



where 



(0) 

Pk 



3fc 



(KV^A + 4a(3A- 1)2) , 



4(3A- 1) 

p;^; = -12fc2(C + 377)^- „;;; (3.6) 



3kV' 



8(3A- 1) 



we find that Eqs. (|3.3[ ) and (|3.4|) can be written in the 
form. 



H 



k SwG 



-Pt, 



a2 3c4 
p, + 3H {pi + p,) = 0, 

where pi = (p, pA, Pfc, Pdr), and 

Pt = ^Pi, pt = ^Pi. 



(3.7) 
(3.8) 

(3.9) 



Eas. (|3.7p and (|3.8p take exactly the forms of those given 
in Einstein's theory of gravity. It is interesting to note 
that the term p^r also appears in the brane world scenar- 
ios BUi. 



IV. THERMODYNAMICS OF THE 
COSMOLOGICAL MODELS IN THE 
HORAVA-LIFSHITZ THEORY 

Thermodynamics in cosmology has been extensively 
studied either in Einstein's theory of gravity [13] or in 
modified theories of gravity, such as brane worlds [2ll |. 
In this section, we shall generalize such studies to the 
HL cosmology. 

The apparent horizon for the FRW model is defined as 

/ ^ g^'f,af,b = 1-(h' + ^) = 0, (4.1) 



from which we find that 



k 



(4.2) 



where f = a{t)r denotes the geometric radius of the two 
spheres t, r = Constants. Following |24j], we define the 
horizon temperature and entropy as. 



Ta 



2'KfA 2tt 



Sa ^ ^ = 



G G{m + M, 



(4.3) 



Introducing the mass-like function M {t, r) [23 and the 
normal vector fc" to the horizon by. 



k" = 5^~rHS^, 



(4.4) 



we find that the energy flow through the horizon is given 
by 



dEA = dE{t,r%^,^ = k'^WaM{t,r%^,Jt 



3/2 ■ 



(4.5) 



g(h + ^ 

On the other hand, from Eq. (|4.3p it can be shown that 



TAdSA 



H H- 



G H 



3/2 ' 



(4.6) 



that is, the first law, TAdSA = dEA, of thermodynamics 
holds on the apparent horizon. It should be noted that 
the above considerations are purely geometric, and were 
not involved with any field equations. Therefore, they 
hold for any metric theories, as mentioned in [20| . 

In the rest of this section, we consider the first law 
of thermodynamics outside the apparent horizon, which 
can be written as [20l |. 

TdSit, V) = Td{sV) = d ip(T)V) + p{T)dV, (4.7) 

where S{T,V) = s{T,V)V denotes the total entropy of 
the system, s{T, V) the entropy density, and T the tem- 
perature. Here we consider that T and V are two inde- 
pendent variables. In other words, we consider a region 
of the universe with a finite radius r. 

Before proceeding further, we would like to note that, 
in [2^ , it was considered the case where the system con- 
sists the whole region inside the apparent horizon f < f a, 
so that the total volume V of the system depends on T. 
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For detail, we refer readers to [2^. In this paper, we shall 
not consider such a possibility. Then, from the condition 
d'^S/dTdV = d^S/dVdT, we find that 



dp _ p + p 
dT " T ■ 

Substituting the above into Eq. ()4.7|) . we find that 

which has the general solution, 

P + P , So 



(4.8) 



(4.9) 



(4.10) 



where sq is a constant and usually set to zero (20| . How- 
ever, here we shall leave this possibility open. The vol- 
ume V is given by 



V = j Vg^d'x 

Air 

= -^a'Vo{r,k). 



47r 



Vl - kr''^ 



-.dr' 



(4.11) 



Inserting Eqs. (|4.8p - (|4.1H) into Eq. (|4.7p . and considering 
the conservation law (|3.8p for pi — p and pi — p, we find 
that 



din I s 



which has the general solution 



4- I si 



3so 



4TTVo{r,k) J ' 



(4.12) 



(4.13) 



where si is another integration constant. Combining 
Eq. (11331) with Eq. (|iTU)) we find that 



p + p = 



(4.14) 



From Eqs. ([X7|) and ([g^ . it can be shown that Eq. ((i7Ti|) 
can be further written as 



„ 2siG fc k IGttG 

H+^^\jm + — = — + ^^pdr 



where 



k = 



2kn'^p.'^K 
"4a(3A- 1)' 



(4.15) 



(4.16) 



Clearly, the first law of thermodynamics holds only when 
condition (|4.15p is satisfied. In other words, it holds only 
for the fluid that satisfies the above condition. This is 
also true in Einstein's theory of gravity, in which it was 
shown that the first law of thermodynamics requires that 
the fluid must consist only three parts [l^ , the cosmolog- 
ical constant, non-relativistic matter, and dark radiation. 
The three parts are coupled each other as [2^ , 



P = PK 



'2'^/PoPK ( — 
\ a 



r+-(?r 



(4.17) 



. CLASSIFICATION OF THE FRW UNIVERSE 
IN THE HORAVA-LIFSHITZ THEORY 

Considering the equation of state given by 

p^wp, (5.1) 



where w is an arbitrary real constant, from Eq. (l3.4p we 
find that, 



P^ Po 



apx 3(i+u.) 



(5.2) 



where po and oq are the integration constants. Since 
Pq represents the energy density when a = ap, we shall 
assume that it is strictly positive pq > 0. Without loss 
of generality, we can always set ao = 1. Then, it can be 
shown that the Friedmann equation (|3.3p can be cast in 
the form [13, [3, 



^a*^ + V{a)^0, 



(5.3) 



where a* = da{t) / d{Hot) , and 

1 / a, 



Via) 



with 



2 V ai+3« 



Po 



+ nk+ n^a" + ^ ) , (5.4) 



3a(3A - l)Hl 
2kC 



a 



dr 



3a(3A - 1)H§ 
4(C + 3r;)fc2 



a(3A - 1)H^ ' a(3A - 1)H^ 

Thus, the acceleration of the universe is given by 



(5.5) 



dV{a) 
da 



^0 



(5.6) 



As mentioned previously, we shall not impose the 
"detailed-balance" conditions, except the condition given 
by Eq. ()2.10p . so that the Friedmann equation ()3.3p has 
the correct coupling sign between the Hubble expansion 
factor and the matter fields. Under such an assumption, 
all the coupling constants appearing in the Lagrangian 
densities (|2.5p are free parameters, so that all the quanti- 
ties defined in Eq. (j5.5p can have any signs, except for 
for which we assume that it is always positively-defined, 
ri„i > 0. When /j = 0, we have Q.^ — ^dr — 0, and 
the corresponding Friedmann equation reduces to that 
of Einstein's theory, studied in detail in [l3|. Therefore, 
in the rest of this paper, we shall assume that k ^ Q. 
Then, it is found convenient to distinguish the three 
cases: f^A = 0, Oa > and f^A < 0. In each of them 
there are seven sub-cases: 



(i) w > 

(iv) w = 
(vi) w — 



(a) w 



(Hi) 



3' 



(v) 



1 

3' 

- 1 < U) < 



1 1 

3<"<3 = 
1 

3' 



-1; (vii) w < —1. 



(5.7) 



In the following we shall consider each of them separately. 
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A. Ha = 



When f^A = 0, Ea. ([0)) reduces to 



2 

from which we find that 



(5. 



V'{a) 



1 



2a2+3» ((1 + 3^") ""^ + 2f^dra'""') , (5.9) 
where V'{a) = dV{a)/da. 



w > 



In this case, from Eq. (|5.9p we can see that when < 
0, the potential has a maximum at 



(5.10) 



where V {a„iax) — 0. When iliir > 0, such a point does 
not exist, and V{a) is a monotonically increasing func- 
tion. Therefore, we shall consider the two cases fldr > 
and ridr < separately. 

Case A. 1.1) ildr < 0: In this case, we have 



Via) 



-oo, 

-nk/2, 



0, 

<■ oo. 



(5.11) 



Thus, depending on the signs of Qk, the potential has 
different behaviors. 

Case A. 1.1. a) fl^r < 0, fife < 0: Then, the potential 
is given by Curve (a) in Fig[Tl from which we can see that 
there exists a point a = a,„ at which we have V{am) = 
0. Thus, in this case if the universe starts to expand 
at the big bang a(0) = 0, it will expand with a < 
until a = am, at which we have d = 0, but we still have 
a = —H^dV{a)/da < 0. So afterwards, the universe 
will start to collapse, until it reaches the point a(ts) — 
again, whereby a big crunch singularity is developed. The 
evolution of the universe is shown schematically in Fig. 

m 

Case A.l.l.b) fl^r < 0, ilk > 0: In this case, there 
exists a critical value fl^^ < for any given flm and fi^, 
which satisfies the conditions 



Via, 
Via, 



max 1 ^dr) — ^i 

max 1 ^dr^ — ^' 



(5.12) 



as shown in Fig. [T] 

When Q.dr < ^dr < 0, the potential is given by Curve 
(b), from which we can see that now Via) = has two 
positive roots, am and Omm, where am > OLmm- In this 
case, the evolution of the universe depends on its initial 
condition. If it starts to expand at the big bang, it will 
expand until a = am and then collapse to a = within 



finite time, whereby a big crunch singularity is developed. 
This is similar to the last case. However, if the universe 
starts to expand at > amin, it will expand forever 
with a positive acceleration d — —H^dVia)/da > 0. It 
is interesting to note that in the latter case a bouncing 
universe is also allowed. For example, if the universe is 
initially collapsing at > amin with d(fi) < 0, then the 
universe will collapse until a = amin ■ Once it reaches the 
point ttmin, where we have d(imin) = and d(tmm) > 0, 
then the universe will turn around, and starts to expand 
acceleratingly without further turning-back, as shown by 
FigEl 

When Qdr = ^dr, the potential is given by the Curve 
(c), from which we can see that now Via) = has two 
degenerate roots, amax = Om = Omm > 0. If it starts to 
expand at the big bang, it will expand until a = amax- 
Since we have V iamax) ^ = V iamax), now a = amax 
represents a stationary point. But, it is not stable, and 
with a small perturbation, it will either collapse to form 
a big crunch singularity at a = or expand forever with 
d > 0. If the universe starts to expand at Oi > amax, 
it will expand forever with a positive acceleration d = 
—HQdVia)/ da > 0, as shown by Fig. [21 If it starts to 
collapse at > amax, the universe will reach the point 
o-max within a finite proper time, and afterwards it will 
stay there forever. However, since now a — amax is not 
a stable point, with a small perturbation, it will either 
start to expand forever or collapse until a big crunch 
singularity is formed at a(is) = 0, as shown by Fig[3) 

When fldr > ^dr^ the potential is always negative, and 
represented by Curve (d). Now the universe will start to 
expand from a big bang singularity at a = forever. But, 
when a < amax it is decelerating, while when a > amax 
it is accelerating. 

Case A. 1.2) fldr > 0: In this case, we have V'ia) < 
0, and as a result, the universe is always decelerating, as 
a(t) is increasing. The potential is given by Figd) 

When flk < 0, the potential is given by Curve (a) in 
Fig. m from which we can see that there exists a point 
am, for which we have Via < am) < 0, where Viam) = 0. 
The universe in this case starts to expand from a big 
bang singularity at a(0) — with d < until it reaches 
its maximal radius am- Afterwards, it starts to collapse 
until aits) = reaches again, whereby a big crunch is 
developed, as shown by Fig. [S] 

When rik > 0, the potential is given by Curve (b) in 
Fig. m from which we can see that the potential is always 
negative, and V'ia) > 0. Thus, in this case the universe 
is expanding from a big bang singularity at a(0) = 
forever. It is always decelerating, as d cx — dV^(a)/da < 0, 
as shown by Fig. [S] 



2. 



In this case, Eq. (l5.8p reduces to 

v«.)^-i(n. + 5 



(5J3) 



6 



Av(a) 




a 



FIG. 1; The potential given by Eq.lUS} for D.A = 0, m > 1/3 

and Qdr < 0. (a) ilk < 0; (b) Qk > 0, D.dr < ^dr', (c) 
t^fc > 0, Qdr = i^dr', and (d) 0,^ > 0, Qdr > i^dr, where 
C = \D.k\/2. 
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FIG. 2; The evolution of the universe with the potential given 
by Eq.JlSJ for Q.K = ^, w > 1/3 and Q.dr < 0. A big bang 
singularity happens whenever a(0) = 0, while a big crunch 
singularity happens whenever a{ts) = 0. 



where fis = ilm + ^dr- 

When fJ^ > and Sl^ < 0, the potential is given by 
Curve (a) in Fig. HI and the corresponding motion of the 
universe is similar to the case of JIa = 0, w > 1/3, fi^r > 
and rjfe < 0. In particular, there exists a maximal 
radius a™. When a < am we have V{a) < and V'{a) > 
0. Thus, the universe in this case starts to expand from 
the big bang at a(0) = with a < until its maximal 
radius Afterwards, it will start to collapse until the 
moment where a{ts) — again, at which a big crunch 
singularity is formed, as shown by the first case in Fig. 

El 

When il^ > and flk > 0, the potential is given by 
Curve (b) in Fig. [H and the motion of the universe is 
similar to the case of JIa = 0, w > 1/3, fi^r > and 



^A = 
w > 1/3 


^dr < 


c 

^dr < ^dr 


^dr - ^dr 
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I/' 

li t 
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'p(t) / 
/II 




\t 
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FIG. 3: The evolution of the universe with the potential given 
by Ea.(|5.8p for f^A = 0, -u; > 1/3 and Qdr < 0. A bouncing 
universe is allowed in either of these two cases. 




FIG. 4: The potential given by Eq.g 
and Sltjr > 0. (a) Qk < 0; and (b) 
|f^fe|/2. 



for ^A = 0, «; > 1/3 
flk > 0, where C = 



flfc > 0, given by the second case in Fig. [5l 

When ils < 0, from Ea. (|5.13[) we can see that the 
potential is always positive for flk < 0, given by Curve 

(a) in Fig. [S] As a result, the motion in this case is 
forbidden. 

When < and flk > 0, there exists a point amim 
for which we have V{a > amim) < 0, as shown by Curve 

(b) in FiglHl Then, the universe in this case will expand 
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FIG. 5: The evolution of the universe with the potential given 
by Ea.(|5.8p for f^A = 0, > 1/3 and Q.dr > 0. A big bang 
singularity happens at a(0) = 0, while a big crunch singularity 
happens at aits) = 0. 
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FIG. 7: The evolution of the universe with the potential given 
by Ea. (|5.13| ) for Q,h = 0, and w = 1/3. A big bang singularity 
happens at a(0) = 0. 




In this case, Eq. (|5.8p reduces to 



dr 



(5.14) 



Thus, depending on the signs of fldr, the potential can 
have different properties. 

Case A. 3.1) ildr < 0: Then, we have 



V{a) 



oo, a = 0, 

-rife/2, a = oo. 



and 



V'{a) = (^iridrl 



(5.15) 



(5.16) 



Therefore, in the present case the potential always has a 
minimum at 



FIG. 6; The potential given by Ea. (l5T3|l for SIa = and 
w = 1/3: (a) Qs < 0, < 0; and (b) fig < 0, fifc > 0, 



where C = \D.k\ /2. 



from a non-zero radius, say, a = > a„i until a = oo 
with d > 0, as shown in Fig. [71 In this case, a bouncing 
universe is also allowed, if it starts to collapse initially 
at Ui > amin with d{ti) < 0. Then, it will collapse until 
it reaches a = amin- Afterwards, it will expand forever 
with d > 0, as shown by the first case in FiglJl 



2 m 



dr\ 



(5.17) 



as shown in Fig. [H 

When ilk < 0, there exists a critical value, fl'^^, given 

by 



1 - 3w 



(5.18) 



so that when \ fldr\ > \^dr\ potential V(a) is always 
positive, as shown by Curve (a), and the motion is for- 
bidden. When Ifldrl = \^dr\ the potential V{a) is always 
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\ \ V ^(c) 
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asX ai \ 
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y(d) 



FIG. 8: The potential given by Ea. ((5TT4)) for ?Lk = 0, ST^dr < 
and -1/3 < ™ < 1/3: (a) < 0, jt^drl > \^dr\; (b) 

fifc < 0, l^drl = (c) Qk < 0, jfidrl < \fldr\; and (d) 

fifc > 0, where C = l^fel /2 and |f2drl is defined by Ea. ((5A5|) . 



positive, except the point a = at which we have 

V{a^in,^dr) = = V'{ajnin,^dr)^ shown by Curve 
(b). This point represents a static universe. In contrast 
to Einstein's static universe in General Relativity [1^, 
this static universe seems stable, as now it corresponds 
to a minimum of the potential. When \fldr\ < l^dA ^^le 
potential V{a) is non-positive only for a € [ai, 02], where 
ai,2 are the two positive roots of V{a) = with 02 > ai. 
Then, the universe is oscillating between the two radii 
oi and 02 without forming any kind of spacetime singu- 
larities. During the period a G [ai,amin) the universe is 
accelerating, while during the period a G {amin,ci2) it is 
decelerating, as shown in Fig. [9l Therefore, in this case 
we have a bouncing cyclic universe. 

When rik > 0, the potential is negative only for a > 03 
where 03 is the real and positive root of V{a) = 0, as 
shown by Curve (d) in FigEl Therefore, in this case 
the universe starts to expand from a non-zero and finite 
radius, say, ai > 03, and shall expand forever. During 
the period 03 > a < amim it is accelerating, while during 
the period a > amim it is decelerating. But the universe 
never stops expanding until a = 00, as shown in Fig|51 A 
bouncing universe is also allowed, if it starts to collapse at 
ai > 03 with a{ti) < 0. Then, the corresponding motion 
is similar to that given by the first case in Fig [31 

Case A. 3. 2) ildr > 0: In this case we have 

na)-(T'/o (5.19) 

^ ^ [ -ilk/2, a = 00, ^ ^ 

and 

V^'(a) = l(^|f^,.| + ii±^^ai-3-) >0, (5.20) 
where equality holds only when a = 00. 



Q.f, = 0. SJj, < 
-1/3 < w < 1/3 


£2 k < 


Qt > 






fid, > 'i'j, 


a(t) 


No Motion 


^ — clmin 




a(t) 




a(l) 
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p(t) 


No Motion 


P =Po 








p(t) 









FIG. 9: The evolution of the universe with the potential given 
by Eq.dSH} for = 0, fidr < and -1/3 <w < 1/3. 

When ilk < 0, the potential has the same form as that 
given by Curve (a) in FigH] for the case f^A = 0, w > 
1/3, ftdr > and ilk < 0. As a result, the motion of the 
universe is also similar to that case. 

When Slfc > 0, the potential is given by Curve (b) in 
Fig. mfor the case 17a = 0, w > 1/3, fldr > and fifc > 0. 
Therefore, in this case the motion of the universe can 
be immediately deduced from there, and given as that 
described in Figl5]for flk > 0. 

It should be noted that although in these two cases 
the motion of the universe has similar characteristics, the 
detail could be different, and when one fits the models 
to observational data, one can get completely different 
conclusions. Since in this paper we do not consider the 
fitting, we shall not distinguish them here. 

4. w^-^ 
In this case, Eq. (|5.8p reduces to 

na) = -l(n^ + ^), (5.21) 

where Ha = Ofe + Qm- 

When ridr > and flk < 0, the potential is given by 
Curve (a) in Figd) Therefore, the motion of the universe 
is similar to that corresponding case. 

When ildr > and flk > 0, the potential is always neg- 
ative and asymptotically approaches to — 12^/2, as shown 
by Curve (b) in Fig. [H Clearly, in this case the universe 
starts to expand from the big bang singularity at a(0) = 
and then expands forever with a < 0. 

When ildr < and Ofc < 0, the potential is always 
positive, as shown by Curve (a) in Fig|Sl So, the motion 
is forbidden. 

When ridr < and Sl^ > 0, the potential is nega- 
tive only when a > amim, as shown by Curve (b) in 
FigEl Therefore, in the present case the universe will 
start to expand from a finite radius, say, ai > am, and 
shall expand forever. Since now we have dV{a)/da < 0, 
we can see that in this case the universe is always accel- 
erating. Note that now no spacetime singularity exists. 
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as we always have a > amim > 0. Certainly, in this case 
a bouncing universe is also allowed. 



-Kw < 



In this case, Eq. (|5.8p can be written as 

1 I ^ ^dr 



Via) = --[nk 



3w -1 



(5.22) 



,3|-uj| + l 



- 2n 



dr 



from which wc have 

^'(a) = -^((3k|-l)a, 

Case A. 5.1) il^r > 0: In this case we have 

F(a) = (-°°' " = 
^ ' L — oo, a = CO, 

and the potential has a maximum at 

_ f^\ndr\_\^' 
-[(3\w\-i))n„J 

as shown in Fig. [TUl 

When rik < 0, there exists a critical value O^,,, 



(5.23) 



(5.24) 



(5.25) 



3|w| + 1 



(5.26) 



so that when fi^r < ^^"^ potential has two positive 
roots, say, ai and 02, where 02 > ai, as shown by Curve 

(a) in FigHni where V{amaxMdr) = V'iamax,^dr) = 0- 
Then, we can see that the motion can have two differ- 
ent kinds, depending on the choice of the initial condi- 
tion of the universe. It can start to expand from the big 
bang at a(0) = until its maximal radius ai and then 
starts to collapse. The collapsing process is exactly the 
time-inverse process of the expansion, and in particular, 
a big crunch singularity is formed at t = ts > where 
a(ts) = 0. Since dv{a)/da > for any given value of 
a € (0,ai], we can see that in this case the universe is 
decelerating, as shown in FigHT] If the universe starts to 
expand at a radius where > 02 , then we can see that 
it will expand forever, and the corresponding acceleration 
is always positive. 

When fldr = ^dr' motion can have two different 
types, too, and the only difference between the last case 
and the current one is that now the two roots ai and 02 
are degenerate and are all equal to Umax, as shown by 
Curve (b) in Fig[TUl As a result, this point represents an 
unstable static point, and any kind of perturbations will 
lead the universe either to collapse or to expand forever, 
as shown in Fig[TlJ When fldr > ^dr^ the potential is 
always negative, as shown by Curve (c) in Fig[Tni and 
the universe will start to expand from the big bang sin- 
gularity at a(0) = forever. Initially the universe is 
decelerating, but once it expands to a^axi it will be ac- 
celerating. 
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FIG. 10: Tlie potential given by Eq. lfOSj) for Qa = 0, Qdr > 
and -1 < ui < -1/3: (a) fldr < ^dr, < 0; (b) Qdr = 
n^ar, > 0; (c) ildr > ildr, ^k < 0; and (d) Qdr > 
0, > 0, where fi^,, is defined by Eq. (f5?^ . 
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FIG. 11: The evolution of the universe given by Ea. (l5.22|l for 

Qa = 0, Qdr > and -1 < w < -1/3. 



When ilk > 0, the potential is always negative, as 
shown by Curve (d) in FigllOl and the universe can start 
to expand from a big bang singularity at a(0) = until 
a = oo. In this case there is no turning point. The 
universe is initially decelerating until a = amax and then 
turns to expand acceleratingly, as shown in Fig. 1111 

Case A. 5. 2) fldr < 0: In this case, from Eq. (|5.22p we 
can see that the potential is always positive for fi^ < 0, 
so that the motion is forbidden. When fife > 0, the 
potential is monotonically decreasing, as shown by Curve 
(b) in FiglSl from which we can see that the potential 
becomes negative when a > amim and a{a > amm) > 0, 
where V{amin) = 0. Thus, in the present case the motion 
of the universe is restricted to a > amin > 0, and no big 
bang or big crunch singularity is developed. The universe 
expands from ai > amin forever until a — oo. Note that 
the spacetime is not singular even at a = cxd, as shown in 
FigEl If initially the universe is in its collapsing phase, 
where Oi > amin and d{ti) < 0, a bouncing universe is 
also allowed. 
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6. w = —1 

In this case, Eq. (|5.8p can be written as 

V{a)=-^ (nk + n,na^ + , (5.27) 

from which we have 

V'{a) = (n,^a^ - ^) . (5.28) 

Then, it can be shown that the potential is given by 
Figim where Curve (a) corresponds to fi^r > 0, fife < 
0, r!™ < Curve (b) to Vldr > 0, f7fc < 0, = f^^; 
Curve (c) to fid,- > 0, ilfe < 0, flm > fl^; Curve (d) to 
^Idr > 0, > 0; and Curve (e) to ildr < 0. 

Comparing Fig. [12] with Fig[TO] we can see that the 
potential has the same shape in each corresponding case 
of (a)-(d), so the motion of the universe in the present 
case can be immediately deduced from there, and the 
corresponding motion of the universe in each case is given 
by Fig. 111! In addition, comparing Curve (e) in Fig. 
[T^ with Curve (b) in Fig. [SI we can see that they are 
similar, except that now the potential goes to — oo as 
a —^ oo. However, this affects only the amplitudes of 
the expansion velocity and acceleration, and the main 
characteristics of the motion are the same in both cases, 
and is given by Fig. [7| However, if the universe chooses 
to collapse initially, a bouncing universe will be created, 
as shown by Fig[3| 

In summary, the motion of the universe in this case is 
given by Figs. [T3| for Curves (a)-(d) in Fig[T2l and by 
the first case in Fig[3|for Curve (e). 

7. w < -1 
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\(e) 









\ ai / 
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FIG. 12: The potential given by Eg. (107)1 for SIa = and 

w = -1: (a) Qar > 0, fi^ < 0, ft„, < ft^^; (b) Qdr > 
0, < 0, n„ = Q"^; (c) fldr >0, nk < 0, n„ > n^; (d) 
fldr > 0, flk > 0; and (e) ftdr < 0. 
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FIG. 13: The motion of the universe with the potential given 
by Eq.((527} for Qa = and w ^ -1. 



In this case, the corresponding fluid is usually called 
phantom, and Eq. (|5.8p can be written as 

Via) = ^l(nk + n,^a^\-'\-' + ^y (5.29) 
from which we have 

V'{a) = (^arr.{3\w\ - l)a^l-l-2 _ (5 

Clearly, when Qdr > the potential has a maximum at 

_ f 2ndr \ ^ 

\(S\w\ - 1)S2,„/ 

Meantime, if flk < 0, then there exists a critical value 
of n^, such that V{amax,^m) = V'{amax,^m) = 0, as 
shown by Curve (b) in Fig[T2| Then, when ri„i < il^ 
the potential will have two positive roots, as shown by 
Curve (a) in Figfl^ while when fl^ > it is always 
negative, and the corresponding curve is that of Curve (c) 



in Fig[T2| When flk > the potential is always negative, 
and is given by Curve (d) in Fig[T2l 

When fldr < 0, the potential is monotonically decreas- 
ing, and is that of Curve (e) in Fig[T21 Therefore, in the 
present case, the motion of the universe is qualitatively 
the same as the corresponding one in the last case, given 
by Fig. [T51and the first case of Fig[3| The only difference 
is that now the matter part is singular as a 00. So, 
now we have a big rip singularity at a = 00, as shown by 
Fig. M 

B. S^A > 

When f^A > 0, from Eq. ([0)l we find that 

Via) = -i(^^. + ^^AaA^ + ^),(5.32) 

= 2a3(i+^) )■ ^'-^^^ 
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FIG. 14: The motion of the universe with the potential given 
by Eq.JS^lJ for Qa ^ and w < -1. 



w > 



When w > 1/3, we first re- write Eas- fOS]) and ([Qa]) 
in the forms, 



in the form, 



3(1 + w)nAa^ + (1 + 3w)f7fe = (3w - 1) 



dr\ 



dr\ 



1 + 3w f2^ 



a3(n-tu) 



(5.38) 
(5.39) 



Clearly, for any given fi^ and flm, the above equations 
always have a solution {a,rt\) — [dmax > 0,0^ > 0). 
Then, we can see that in the current case we also have 
three different sub-cases according to whether I^a < fi^, 
r^A = ^% or r^A > f^Ai which the potential is given, 
respectively, by Curves (a), (b) and (c) in Fig. [TOl and 
the motion of the universe is given in each case by the 
corresponding one given in Fig. Illl and Fig|31 



V{a) 



V'{a) = -airiA- 



1 



a3(i+«)) 

X I i^^r^r!™ + r!<i,a3— 1 ) ) . (5.35) 



To study the solution further, we consider the cases 
fldr > and ildr < separately. 

Case B.1.1) fldr > 0: In this case, from Eq. (|5.34p 
we find that 



Via) { 



— CX3, a = 0, 
—oo, a = cxD, 



(5.36) 



while from Ea. (|5.35p we can see that the potential has a 
maximum at amax, where V'{amax) — 0. 

Case B.l.la) Qdr > 0, fife < 0: If fife < 0, from 
Ea. (|5.34p we find that we can further have V{amax) ~ 0, 
that is, for any given flk < and Cldr > 0, there always 
exists a positive solution, {amaxT^%), of the equations. 



In this case, we have, 



Via) = -- [^k + ^Ao' 



a2 



(5.40) 
(5.41) 



where JIa = ^^r?j + ^dr- 

Case B.2.1) i^A > 0: In this case we find that the 
potential always has a maximum at 



riA 



1/4 



for which we have 



(5.42) 



(5.43) 



Thus, if rifc < 0, a critical point Q\ exists, so that 
V ia,nax, ^a) = V ia^nax, = 0, whcrc 



4r!A 



(5.44) 



1 + 3u; 



3(l+tu) 
amax 



ildrO, 



max J 



3w-l 
max 



(5.37) 



When r^A < i^Ai potential will take the form of 
Curve (a) in Fig. [TOl while when fl\ = fl'jy or fl^ > f2^, it 
will take, respectively, the form of Curves (b) and (c), so 
the motion of the universe is similar to the corresponding 
cases given in Fig. 1111 

Case B.l.lb) ildr > 0, > 0: In this case, it can 
be shown that the potential is given by Curve (d) in Fig. 
1101 so the motion is given by the corresponding case in 
Fig. M 

Case B.1.2) fldr < 0: In this case, it can be shown 
that the equations Via) — and V'ia) = can be cast 



Therefore, in this case, there are three sub-cases accord- 
ing to whether f^A < ^^A' ~ '^^ > ^A' 
which the potential is given, respectively, by Curves (a), 
(b) and (c) in Fig. [TOl and the motion of the universe 
is given in each case by the corresponding one given in 
Figs, [n] and [3] 

When rik > 0, V'ia) is always negative and the poten- 
tial is given by Curve (d) in Fig. [TOl and the motion of 
the universe is given by the corresponding case given in 

Fig. m 

Case B.2.2) J^a = 0: In this case we have 



Via) 



i {nk + nAa^) 



(5.45) 



and the potential is given by Fig. [TOl from which we 
can see that when flk < 0, Via) is negative only when 
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V(a) 




FIG. 15: The potential given by Ea. (fO0|l for lu = 1/3. (a) 

J^A > 0, Oa = 0, fife < 0; (b) > 0, ^a = 0, fife > 0; and 
(c) > 0, r^A < 0, where C = |r2fe| /2. 



a > arm where V{am) — 0, as shown by Curve (a) in 
FigHH Thus, the motion in this case is restricted to 
a > am- Since dV{a)/da < 0, we can see that now 
the universe is always accelerating, as shown in Fig ll6l 
When flk > 0, the potential is always negative [as shown 
by Curve (b) in Fig. 1151 and the universe can start to 
expand from the big bang at a(0) = until a = oo. 
Case B.2.3) 51a < 0: In this case we have 



Via) 
Via) 



1 



2 

-a ( il\ 



A 



a^ 
< 0. 



(5.46) 



Therefore, in this case the potential is a monotonically 
decreasing function, as shown by Curve (c) in Fig |15l 
Then, we can see that the motion is restricted to the 
region a > a„j. The universe is always accelerating, and 
no singularity exists in the present shown in Fig. 

1161 A bouncing universe is also allowed in the present 
case, and the corresponding motion is described by the 
first case in FiglSl 



3. -1/3 <w <l 

In this case we find that the potential and its first 
derivative can be written as 



Via) 



1 



(5.47) 
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FIG. 16: The motion of the universe with the potential given 
by Eq.fOni for ^a > and lo = 1/3. 



Via) 



-a I Q\ 

X [rtdr- 



1 

l + 3w 



^rr}.a 



l-3w 



(5.48) 



To study this case further, we need to consider the cases 
fldr > and ildr < separately. 

Case B.3.1) fldr > 0: In this case from the above 
equations we find that 



Via) = { 



-oo, a — 0, 
-oo, a = oo, 



(5.49) 



and that the conditions Via) — V'ia) — can be written 
as 



l + 3w 



SIa = — 7 ^dr 

a^ \ 



l + 3w 



1 - 3w 
2a2 

l-3w 



rid,, (5.50) 
(5.51) 



Clearly, for any given Qdr > and flk , the above equation 
always have a positive solution (a, Q\) = iamax, ^a)^ ^'^ 
that when Q\ < SIa, the potential is given by Curve (a) 
in Figlini when Ha = il'j^, it is given by Curve (b), and 
when r^A > f^A' given by Curve (c). Then, we can see 
that the motion in this case is given by the corresponding 
case given in FigfTTl 

Case B.3.2) ildr < 0: In this case from Eqs. (l07|) 
and (|5.48p wc find that 



Via) = { 



oo, a — 0, 
— oo, a — oo. 



(5.52) 



and 
Via) 



.1 (^r!,a4_i±^0 



-^Fia) 



(5.53) 
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FIG. 17: The function F{a) defined by Eg. dKSSll for > 
0, fidr < and -1/3 < < 1/3. (a) \^ir\ > \^dr\; (b) 
\Qdr\ = inS.l; and (c) in^.l < 



From the above expression we can see that there exists 
a critical value of fJ^^, so that when \fldr\ > \^dr\' 
function F{a) is always positive, and V'{a) = has no 
real positive root, as shown by Curve (a) in FiglTTl When 
\fldr\ = \^dr\^ equation F{a) = has only one real 
positive root, as shown by Curve (b), and when \ fldr\ < 
\^dr\^ ^i'^) — has two real positive roots, as shown by 
Curve (c) in FigfTTl 

On the other hand, the conditions V{a) = V'{a) = 
can be written as 



l + 3w^ , 3(w + l)_^ 2 



Ha 



-ilk 



(1 + 3u;)f7,; 



1 - 3w 
2a? 



0,(5.54) 
(5.55) 



From Eq. (|5.54p we can see that it has solution only when 
rife < 0. Thus, in the following we need to consider the 
cases ilfe < and ilfc > separately. 

Case B.3.2a) ilfc > 0: In this case from the above 
we can see that when \ildr\ > \^dr\' potential is a 
monotonically decreasing function, and given by Curve 



riS I, the potential has 



(d) in Fifdll When {ildrl 
only one minimum, and there are three different cases, 
as shown by Curves (a)-(c). In each of these four cases 
(a)-(d), we can see that the potential is negative only for 
a > ttm where a™ is the positive root of ^(a) = 0, as 
shown by Fig. [THl Therefore, the motion of the universe 
is restricted to a > am- If the universe starts to expand 
from ai > a^ it will expand forever. Since dV(a)/da < 
for a > arm we can see that the universe in each of these 
cases is accelerating. If it is collapsing initially, then it 
will reach its minimal radius am within a finite proper 
time. Afterwards, it will start to expand, and a bouncing 
universe is produced, as described by the first case in 
FigEl 

When \Vldr\ < \^dr\^ ^^'^ potential has one mini- 
mum and one maximum, as shown by Curve (e) in 
Fig. [THl This is a very interesting case. As Cases 
(a)-(d), the motion of the universe is also restricted to 




FIG. 18: The potential given by Eg. dOT)) for -1/3 < to < 
1/3 and Qdr < 0, Qk > 0. (a) l^drl = \iidr\ , V'(a,„i„) > 

0; (b) l^drl = 1^5,1, V{armn) = 0; (C) l^drl = 

in^,,!, V{a^in) < 0; (d) l^drl > |!^drl, V{amin) > 0; and 

(e) \ndr\<mr\- 



the region a > am- But, it is fundamentally differ- 
ent from these cases: The universe is accelerating for 
a e [am, cLmin) and a G {amax, oo), and decelerat- 
ing for a S {amim CLmax)- Therefore, it can describe 
the evolution of our universe without a big bang singu- 
larity. In particular, if the universe chooses to collapse 
first at ai, where amin < ai < amax, we can see that a 
scale-invariant perturbation can be produced during this 
matter-dominated period [131 , and the universe will expe- 
rience a bouncing once it collapses at am, whereby a big 
bang singularity is avoided. Once it turns to expand, it 
will first expand acceleratingly until a — amin- Clearly, 
if this expansion is large enough, the horizon problem 
can be solved. Afterwards, the universe will experience 
a decelerating period until a = amax- Once this point 
reaches, it will expand with a positive acceleration, which 
may be identified with the late cosmic acceleration. 

Case B.3.2b) ilk < 0: Now it can be shown that in 
the sub-cases \ildr\ > \^dr\ ^"^^ \il 



1dr\ = \^dr\ poten- 
tial is quite similar to the corresponding cases given in 
FigdU that is, it is given by Curve (d) for Ifldr] > \^dr\ 
and Curves (a)-(c) for Ifldr] — \^dr\- However, when 
\ildr\ < \^dr\: Conditions V{a) V'{a) have pos- 
itive root for a, and now we have five different cases, as 
shown in FigdU In the case described by Curve (a) the 
motion of the universe is restricted to the region a > am, 
and it can represent an expanding or a bouncing uni- 
verse, depending on the initial velocity of the universe. 
In the case described by Curve (b), the motion of the 
universe is similar to the last case, except that now a 
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FIG. 19: The potential given by Eg. ((07)1 for Sldr < 0, flk < 
0, i^drl > Iffa^l and -1/3 < w < 1/3. (a) l/(a„,„) > 
0, V(a,„„a;) > 0; (b) T/(o™„) = 0, V{a„,ax) > 0; (c) 

V(a„i„) < 0, V{amax) > 0; (d) V{amin) < 0, l^(amaa:) = 0; 

and (e) V'(ami„) < 0, V{amax) < 0. 



stationary universe also exists at a = dmin- In the case 
described by Curve (c), the motion of the universe for 
a > am is similar to the last two cases, but now a bounc- 
ing cyclic universe exists for a € [am, aj- For the case 
described by Curve (d), if the universe starts to expand 
at ai > am, it will expand forever with a > 0. If it col- 
lapse from a point > am, once it reaches a — amax it 
will stay there. However, since it is a non-stable point, 
with a small perturbation, the universe either continu- 
ously collapses until a = am, or starts to expand forever. 
If it continuously collapse, when it reaches am, it will 
start to expand until amax- The following motion can 
either continuously expand or collapse. In tis case, the 
universe can also move between a,„ and amax, as shown 
in Fig[T9l In the case described by Curve (e), the motion 
is the same as the corresponding case in FigfTSl where a 
bouncing cyclic universe is produced. 



4. w = -L 

In this case we find that the potential and its first 
derivative can be written as 

V{a) = -\{vLi^ + ^li^a^ + (5.56) 



V'{a) 



A / 4 



(5.57) 



To study this case further, we need to consider the cases 
Q.dr > and fldr < separately. 



Case B.4.1) fldr > 0: In this case from Ea. (|5.56p we 

find that y(a) — > — oo as a — > and a —^ oo, while from 
Eq. (|5.57p we can see that the potential has a maximum 
at amax — (f^dr/f^A)^^^, for which we have 



V{a„ 



1 



(5.58) 



Therefore, if JIa < 0, there exists a critical value of 
il^ = ri^/(4i7dr, for which we have V {amax,^']\) — 
V {amax,^\) = 0, as shown by Curve (b) in FigfTOl 
When a. A < the potential is given by Curve (a), and 
when < Q\it is given by Curve (c). When JIa > 0, 
the potential is always negative, and described by Curve 
(d) in Fig[TOl Thus, in these sub-cases the motion of the 
universe is described by the corresponding cases given in 
FigdU 

Case B.4.2) ri^r < 0: In this case, we have 



Via) = { 



oo, 



a 
a 



0, 
oo, 



(5.59) 



and V'{a) is always negative. Then, the potential is given 
by Curve (e) in Fig[T2l and the corresponding motion of 
the universe is described by the corresponding case given 
in Fig[T3] and the first case of Fig [3] The only difference 
is that now the matter density diverges at the big bang 
singularity a(0) = 0. 



5. -Kw < -1/3 

In this case we find that the potential and its first 
derivative can be written as 



Via) 



Via) = 



+a3l'"l-i(an + f^Aa'^^~l"l))), (5.60) 



(5.61) 



To study this case further, we need to consider the cases 
Q.dr > and fldr < separately. 

Case B.5.1) ftdr > 0: In this case from the above 
equations we find that 



Via) = ( " " ' (5.62) 

and that the conditions Via) — V'ia) ~ can be written 



as 



1 / 

flk = 2 [ + a' 

( 



a^ 

X (ri„ 



a 



3|i«| + l 
3(1-1-1)^^^ 



(5.63) 



f^A = ^ Udr - . (5.64) 
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Clearly, they have solutions for positive a and only 
when < 0. Then, we have three different sub-cases, 
as given by Curves (a), (b) and (c) in Fig. [TOl which 
correspond to, respectively, fl^ < O^, n\ = and 
ft A > S^A, where fi^ is the solution of Eas. (|5.63[) and 
()5.64p . When fife > the potential is always negative 
and is given by Curve (d) in Fig. [TOl so the motion of 
the universe is described by the corresponding cases in 
FigHH 

Case B.5.2) ildr < 0: In this case we find that 

y(a) = (°°' (5.65) 

and V'{a) < 0. Therefore, now the potential is mono- 
tonically decreasing, and is given by Curve (e) in FigfT2l 
Then, the motion of the universe in this case is given 
by the last case described in Fig[T31 except that now we 
still have p(0) = oo and p{oo) = 0. Certainly, a bouncing 
universe is also allowed. 



6. w = —I 

In this case, we find that 

V{a) = (nk + nmAo" + ^) , (5.66) 

V\a) ^ -a (^f7„,A - ^) , (5.67) 

where rim, A = + f^A > 0. Comparing the above 
equations with Eqs. (|5.27p and (|5.28p in the case J^a — 0, 
we can see that by exchanging f2„j with f2m,A, we can 
get one set of the equations from the other. Therefore, 
the motion of the universe in the present case can be 
immediately deduced from the corresponding ones given 
in Figs. [11 [13 and [3 



and that the conditions V{a) = V'{a) = can be written 
as 

f2fc = -a2(l]A + a'l'"l"'Om) (5.71) 

nA + ^ + ^J4-^"'""""""'- ^^-^2) 

Clearly, they have solutions for positive a and Oa only 
when rjfc < 0. When Vlk > 0, the potential is strictly 
negative, and is given by Curve (d) in Fig[T21 When 
fife < 0, Eqs. (|5.7ip and (|5.72p have a unique solution 
iamax,^%), SO that the potential is given, respectively, 
by Curves (a), (b) and (c) in Fig. [TOJ which correspond 
to JIa < f^A' = f^A' > f^A- Therefore, the 

motion of the universe can be deduced from there and 
can be shown that it is given by FigfMl 

Case B.7.2) ildr < 0: In this case we find that 

y(a) = (°°' (5.73) 
^ ' L — oo, a = oo, ^ ' 

and V'{a) < 0. Therefore, now the potential is mono- 
tonically decreasing, and is given by Curve (e) in FigfTOl 
Then, the motion of the universe in this case is given 
by the last case described in Fig[T3] and the first case in 
FigEl The only difference is that now the spacetime is 
singular at a = oo. 

C. Qa <0 

When JIa < 0, from Eq. ([0)) we find that 
Via) . -l{n,-\nA\a^ + ^ + -^yi5.7A) 

X (^l±^n„^ + a^--'nd,^y (5.75) 



7. w <-l 

1. w>^ 

In this case we find that the potential and its first 
derivative can be written as In this case, we have 



+aMrjA + r2ma^^l"'l"^)]] , (5.68) When n^r > 0, from Eq. (|5.75p we can see that 

^ dV{a)/da is always positive, and the potential is mono- 

y/j-^^ ^ — a ( ilA — ^'^^ tonically increasing, as shown by Curve (a) in Fig. [201 

\ a* Thus, the motion now is restricted to a < a„i. 

3|^| _ 1 I I \ When Qdr < 0, the equations V{a) = and V'{a) = 

+ f^^a'^l-^l-M . (5.69) 

can be written in the forms, 

Case B.7.1) ildr > 0: In this case we find that "fc = ; a H — —^r—, (5.77) 

1 + 3w (l + 3w)a^ 

na) = r°°' (5.70) ^5 78) 

^ ' L — oo, a = oo, ^ ' ' ' 2a-^(i+'") a'* 
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(c) (d) 




FIG. 20; The potential given by Eg. lfOill for w > 1/3. (a) 

fldr > 0; (b) flar < 0, fife > 0, |!^a| > (c) S^dr < 

0, > 0, |!^a| = l^^v!; (d) 0*- < 0, fifc > 0, I^a] < \ni\- 
and (e) fldr < 0, Slfc < 0. 



which have solutions for a > and CIa < only when 
fife > 0. In the latter case, there exist a critical value 
|r2^|, for which, when jr^Al > the potential is given 



by Curve (b), when = I^^aI' by Curve (c), and when 
\nA\ < \ni\, by Curve (d), in FigEOl When flk < 0, 
Eas. (|5.77p and (|5.78p have not positive solutions for a and 
and the potential is given by Curve (e) in FigEOl 
In all the above cases, we can see that the motion is 
always restricted to a < where am is the unique posi- 
tive solution of V{a) = 0, as shown in Figl^Ol Therefore, 
now the universe starts to expand from the big bang at 
a(0) = to its maximal radius a,„, and then starts to 
collapse until a big crunch singularity is formed at the 
moment tg where a{ts) = 0. Note that for the poten- 
tial given by Curves (a) and (b), the universe is always 
decelerating. For the one given by Curve (c), the point 
a = am is a stationary point, as now we have a = = a 
at this point. However, as shown there, this point is not 
stable, and with a small perturbation, the universe will 
collapse towards a = and finally a big crunch singular- 
ity is formed there. For the potential given by Curves 



(d) and (e) , there exist a point a 



ar, 



< am, at which 



a = —HQdV{a)/da = 0, as shown by these curves. How- 
ever, the universe will immediately turns to decelerate. 



2. = ! 



In this case, we have 



Via) = -1 (n,-\nA\a^ + ^^ , (5.79) 
Via) = a (\Qa\ + ^) . (5.80) 



Case C.2.1) < 0: In this case, we have 

y(a) = (°°' " = (5.81) 
^ ' I oo, a = CO. ^ ' 



When rifc < 0, the potential is always positive, as shown 
by Curve (a) Fig. [511 and the motion is forbidden. 
When Jlfc > 0, Via) = = V'ia) have the solution. 



1/2 



(5.82) 



If jr^Al > I^^aI' the potential is positive, as shown by 
Curve (b) in FigETl and the motion is forbidden. If 
I^^aI = I^^aI' the potential is also positive, except for 
the point a = ami„, at which we have Viamm) = 0, as 
shown by Curve (c). At this point we have ct = = a, 
and it represents a stable stationary point, similar to the 
case described by Curve (b) in FigO If I^IaI < I^^aL 
the potential is negative only for the range a € (01,02), 
as shown by Curve (d) in Fig |211 which is also similar 
to Curve (c) in FigjSl and the motion of the universe is 
described by the corresponding case in FiglQ] 
Case C.2.2) fls > 0: In this case, we have 



Via) = - { 



00, 

CO, 



0, 

CO, 



(5.83) 



and V'ia) > 0, that is, now the potential is a mono- 
tonically increasing function, as shown by Curve (e) in 
Figim which is similar to Curve (a) in Figl^Ol There- 
fore, the motion of the universe in this case is similar to 
that case. 



3. -1/3 <w <l/3 
In this case we find that 

Via) = -i(f7fe-|f7A|a' 

1 



V'ia) = a |f7A|+ 4 

X ( Odr + 



a 

l + 3w 



Case C.3.1) Qs > 0: In this case, we have 
Via) = -{: 



CO, a = 0, 
CXD, a — OO, 



(5.84) 



(5.85) 



(5.86) 



and V'ia) > 0, so that the potential now is a mono- 
tonically increasing function, as shown by Curve (e) in 
Figim Therefore, the motion of the universe in this case 
is similar to that one. 
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Thus, when > we have V'{a) > and 



FIG. 21: The potential given by Eg. ((579]) for w = 1/3. (a) 

flk <0, fls < 0; (b) fife >0, ns< 0, |J7a| > I^aI; (c) flk > 

0, ^5 < 0, i^aI = (d) >o, ns < o, |J7a| < 

and (e) > 0. 



Case C.3.2) fls < 0: In this case, we have 

l/(a) = (°°' " = (5.87) 



and V{a) = = V'{a) can be written as 



ii~3w)\ndr\ 

(1 + 3u;)a2 
(1 + 3w)il„i 



Hl + w)\nA\^,^ (5.; 



(1 + 3m;) 



(5.89) 



Clearly, for any given flk and firm the above equations 
always have a positive solution {amin, ^^a)' that when 
ir^Al > \^%\, the potential is always positive, and the 
motion is forbidden, as shown by Curve (b) in FigEI 
When IJIaI = I^aIj the only possible motion is a = amm, 
at which we have a = = a, as shown by Curve (c). 
When IJIaI < I^aI' the potential is negative in the range 
a € (oi, a2), as shown by Curve (d) in FigdlJ Therefore, 
the motion of the universe in the current case can be 
deduced from the corresponding ones given in Fig |21l 



4- w - 

In this case, we find that 



-1/3 



Via) = -Un^^\nA\a' + ^ 



V'ia) = ^ f a* 



n 



dr 



(5.90) 
(5.91) 



na) = {. 



a 
a 



0, 
oo. 



(5.92) 



that is, the potential now is a monotonically increasing 
function, given by Curve (e) in Figl2H and the motion of 
the universe is restricted to a < am- The universe starts 
to expand at the big bang a(0) = until its maximal 
radius am, and then starts to collapse until a big crunch 
singularity is formed at the moment ts where is given 
by a{ts) = 0. 

When ridr < 0, we have 



cxD, a = 0, 



and the potential has a minimum at 



dr 



nA 



1/4 



at which we have 



V [amin) = i^A — 2 Ir^dr^^Al) ■ 



(5.93) 



(5.94) 



(5.95) 



Clearly, if Oa < 0, the potential is always positive, and 
the motion is forbidden. When f^A > 0, there exist a crit- 
ical value f2A; such that the potential is always positive 
for jr^Al > I^^aI' the motion is forbidden, as shown 
by Curve (b) in FigEll When IQa] = \^a\: the only 
possibility of the motion is stationary point a = amin, as 
shown by Curve (c) in FigET] When |17a| < the 
potential is negative in the range a G (ai,a2), as shown 
by Curve (d) in FigEll 



5. -l<M<-l/3 

In this case the potential and its derivative can be writ- 
ten as 



Via) 



1 



O I I „3|«i|-l 

"fc H T' + « 



(a^-|0A|a3(i-H))) 



Via) = a[\nA\ + — - ^^3(,_|,|) 
Case C.5.1) Vldr > 0: In this case, we have 

F(a) = f-°°' " = '^' 
^ ' L oo, a = oo, 

and Via) ^ V'ia) = yield 

l]fe + anfl'l'"!"' + % = \nA\a^, 
a'' 

I , ^d7- (3|w| - l)f7,„ 



2a3(i-K"l) 



(5.96) 
. (5.97) 

(5.98) 

(5.99) 
(5.100) 
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Clearly, for any given Qk and Qm, there always exists 
a solution (a, |r2A|) — (omm > 0, |r2^|), so that when 
ir^Al > the potential is given by Curve (b) in Fig|20l 
when jr^Al = I^aI' given by Curve (c), and when 
|S1a| < I^^aI' given by Curve (d). Then, the mo- 
tion of the universe in this case can be deduced from the 
corresponding cases given in Fig l20l 

Case C.5.1) fldr < 0: In this case, we have 

y(a) = (°°' (5.101) 
^ ' I oo, a ~ oo, ^ ' 

and V{a) ~ V^'(a) — yield the same equations (|5.99[) 
and ()5.100p but now with < 0, which show that for 
any given f2fe, fl^r < and f2„i, they also have a unique 
solution, (a, |17a|) = (flmin > 0, |r2^|),but now when 
I^IaI > I^^aIi potential is always positive, as shown 
by given by Curve (b) in FigET] When jf^Al = I^^aL i* 
is given by Curve (c), and when jr^Al < I^^aI' given 
by Curve (d) in Figl2T] Then, the motion of the universe 
in this case can be deduced from the corresponding cases 
given there. 



6. w = —I 

In this case we have 

V{a) = -^(nk + nm,Aa^ + ^y (5.102) 

where flm,A is defined in Eqs. (|5.66p and (|5.67p . but now 
can be positive, zero or negative. We also have the pos- 
sibilities, f2(jr > and fl^r < 0. In the following we 
consider each of these cases separately. 

Case C.6.1) rim, A > 0, > 0: In this case, from 
Ea. (|5.102p we find that that 

V{a)^\'°°' (5.104) 
^ ' L —oo, a — oo, ^ ' 

and V'{a) — has the solution, amin = i^dr/i^m,Ay^^ , 
at which 

V {a^ax) = -\ (^k + 2 {\n^,Andr\f'^) ■ (5.105) 

Thus, if rifc > 0, we have V {amax) < 0, and the po- 
tential is strictly negative, as shown by Curve (d) in 
FigdU If Vlk < 0, we have V {a^ax,^''m,A) = 0' where 
A = ^fc/(4f^dr, as showu by Curve (b) in FigfT^ 
When nrn,A > ^mA' potential is given by Curve 
(c) there, for which we have V{a) < for any given a. 
When rim, A < ^mA' potential is given by Curve 
(a) in FigHH from which we can see that V{a) < is 
only possible when a € (01,02) where oi 2 are two pos- 
itive roots of V{a) = 0, as shown there. Therefore, in 



this case the motion of the universe can be deduced from 
FigdU 

Case C.6.2) fim^A > 0, ildr < 0: In this case, from 
Eqs. (|5.102p and (|5.l63p we find that that 

F(a) = |+°°' " = °' (5.106) 

^ ' L —00, a = 00, ^ ' 

and V'{a) < 0, that is, now the potential is monotonically 
decreasing, and the motion of the universe is possible only 
in the range a > a„iin, as shown by Curve (e) in Fig |12l 
Case C.6.3) rim,A — 0, ^Idr > 0: In this case, from 
Eqs. (|5.102p and (|5.l63p we find that that 

^w-{:S/2, (^.lor) 

and V'{a) > 0. Thus, now the potential becomes mono- 
tonically increasing. When il^ < 0, the potential is given 
by Curve (a) in Figdl from which we can see that the 
motion of the universe now is restricted to a < am, 
where V{am) = 0. When ilk > 0, the potential is 
given by Curve (b) in FiglH in which the potential is 
always negative, and the universe starts to expand from 
a big bang until it reaches a = 00 with an infinite 
proper time. The motion is always decelerating as now 
a = -HldV{a)/da < 0. 

Case C.6.4) flm,A — 0, i^dr < 0: In this case, from 
Eqs. (|5.102p and (|5.l63p we find that that 

^(«)-{"i7./2, (^-los) 

and V'{a) < 0. Then, it can be shown that when il^ < 0, 
the potential is always positive, and given by Curve (a) 
in FiglG] Therefore, in this case the motion is forbidden. 
When fifc > 0, the potential is given by Curve (b) in 
FigEl in which the potential is negative only when a > 
o-mim the universe is always accelerating, starting from a 
non-singular point > amin- The universe is also free 
from singularity at a = 00. 

Case C.6.5) ilm,A < 0, fi^r > 0: In this case, we 
find that that 

^(«) = {f7j2, aZt, (5-109) 

and V'{a) > 0. It can be shown that for any given fi^, 
the potential is given by Curve (a) in Figl20l Then, the 
motion of the universe can be deduced from there. 

Case C.6.6) ilm.A < 0, ildr < 0: In this case, we 
find that that 

V(a) = (o'/„ (5.110) 
^ ^ [ S2fe/2, a = 00, ^ ' 

and V'{a) — has the solution amin = \^dr/^m,A\^^'^ i at 
which we have 

V {am^n) = -\ (r!fc - 2 {\^mA^dr\i''-^ ■ (5.111) 
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Thus, if flk < 0, we have V {umin) > 0, and the poten- 
tial is strictly positive, as shown by Curve (a) in Figl^TJ 
So, the motion now is forbidden. If Qk > 0, we have 
y(a™„,f7f„ ^) = 0> where |f^™,A| = f7^/(4 jf^rf^l, as 
shown by Curve (c) in FiglJTJ Then, the only possible 
motion is that the universe is static and stays at the 
point a = amim at which we have a = a = 0. So, 
it represents a stable point. When |ri„i,A| > |f^m^A|' 
the potential is always positive, as shown by Curve (b) 
in Fig[2Tl so the motion in this subs-case is forbidden. 
When |ilm,A| < I^toaIj the potential is given by Curve 
(d) in Fig|2n from which we can see that the motion is 
restricted to the region a g [ai , 02] , and the universe is 
oscillating between this two turning points, and no space- 
time singularities is formed during the whole process. 

7. w <-l 

When w < -I, Eqs. ((??7i)l and can be written 

as 

-a2(|r!A|-a3(l-l-i)r!„^)), (5.112) 
V'{a) ^ a{\n^\ + ^ 

_(^H_l)a^^3(|H-i)^ . (5.113) 

Case C.7.1) Vldr > 0: In this case we have 

V{a) = \'°°' " = (5.114) 

It can also be shown that the equations V{a) = = V'{a) 
have a unique solution {am,ax,^\)- The potential in the 
cases, l^iAl > I^^aI = l^^^l, and |r!A| < is 

given, respectively, by Curves (a), (b), and (c) in Fig. 
1121 and the corresponding motion of the universe can be 
easily deduced from there, as described by FigfT4l 
Case C.7.2) fldr < 0: In this case we have 

y(a) = (°°' (5.115) 
^ ' I —00, a = CXI, ^ ' 

and the equations V{a) = = V'{a) have a unique solu- 
tion (amaxT^A)- The potential in the cases, |r2A| > I^^aI' 
ir^Al = I^aI' s-iid I^^aI < I^^aI' is given, respectively, by 
Curves (a), (b), and (c) in Fig. [THl and the corresponding 
motion of the universe can be easily deduced from there. 

VI. CONCLUSIONS 

In this paper, we have studied the thermodynamics 
of cosmological models in the Horava-Lifshitz theory of 
gravity, and found that the first law of thermodynamics 



holds only for the perfect fluid that satisfies the condition 
(SIB. 

In the Horava-Lifshitz theory of gravity, the 
Friedmann-like equations are given by Eas. (|3.3p and 
(13. 4p . We have studied systematically these equations 
coupled with a perfect fluid with the equation of state 
p = wp, where p and p are the pressure and energy den- 
sity of the fluid, and w is a constant. In this case, the cor- 
responding cosmological models contain four free param- 
eters, r2„i, flk, ^dr and fl^, where both f2fe and ^dr are 
proportional to the curvature of the three-dimensional 
space [cf. Eq. (|5.5p ]. and arc all zero when the curvature 
vanishes. Then, the models reduces to those given in 
Einstein's theory of gravity, which has been systemati- 
cally studied and classified recently in [l3|. Therefore, in 
this paper we have assumed that flk^dr 0, but kept 
their signs arbitrary. The term f^A is related to the cos- 
mological constant, and we have studied all the three 
possibilities, f^A = 0, f^A > 0, and f^A < 0, although we 
have assumed that il™ > 0, as the latter represents the 
matter. Then, depending on particular values of the four 
free parameters, we have divided the models into various 
cases. In each case the main properties of the evolution 
have been studied in detail, including the periods of de- 
celeration and/or acceleration, and the existence of big 
bang, big crunch, and big rip singularities. 

As first noticed in 10], models that represent a bounc- 
ing universe can be constructed, due to the presence of 
the dark radiation term. However, it should be noted 
that the condition fJ^r < for the existence of such mod- 
els is only a necessary condition, but not sufficient. For 
example, in Case c) where ^a < 0, no such models exist 
for all the sub-cases with w > —1. In addition, in order 
to have a bouncing universe, it is also important that w 
has to be less than 1/3, so no source of matter will red- 
shift faster than that of the dark radiation. However, as 
argued in [4], the radiation in the UV regime scales as 
a~^. Thus, in order to obtain a bouncing universe, one 
mig ht need to consider the more general case studied in 
[ISj . in which a term that proportional to a^^ exists. So, 
if this term dominates the radiation, a bouncing universe 
still exists. 

It should be also noted that many of these models may 
not be consistent with current observational data sets 
As a matter of fact, if we compare them with the 
ACDM model, we find that 

VACDM{a) = (^^k + f^Afl^ + 

^ r-00, a = 0, 

L —00, a = 00, ^ ' 

which is schematically given by one of the curves given 
in Figlinior FigfT2l except for the case of Curve (e) given 
there. From these curves we can see that there exists 
a maximal point a = amaxi for which when a < amax 
the universe is decelerating, and when a > amax it is 
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accelerating. From Eq. (|6.ip we find that 

Since our universe now is accelerating, and recall that 
we have set the current radius ao of our universe to one, 
from the above expression we can see that the current 
acceleration of the universe requires that < 2n\. In 
fact, numerical fitting of the model shows that flm — 0.27 
and — 0.73 (28| . for which we have Umax — 0.57. 
However, we do hope that such a classification is useful 
for the future studies of cosmology in the Horava-Lifshitz 
theory, as such studies are just starting 

Finally, we note that in brane world scenarios , dark 



radiation term also appears. In the later times of the 
universe, the quadratic terms of the energy density can 
be neglected, and the corresponding Friedmann equation 
will reduce to the ones studied in this paper. So, our 
classification presented here is also applicable to these 
models, too. 
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